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Magnetic fields created in the noncentral heavy-ion collision are studied within a mi- 
croscopic transport model, namely the Ultrarelativistic Quantum Molecular Dynamics 
model (UrQMD). Simulations were carried out for different impact parameters within 
the SPS energy range (Ei a b = 10 — 158A GeV) and for highest energies accessible for 
RHIC. We show that the magnetic field emerging in heavy-ion collisions has the magni- 
tude of the order of eB y ~ 10 _1 • m\ for the SPS energy range and eB y ~ for the 
RHIC energies. The estimated value of the magnetic field strength for the LHC energy 
15 -m%. 
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1. Introduction 

Heavy-ion collisions are intensively investigated both experimentally and theoret- 
ically to reveal information about properties of the nuclear matter under extreme 
conditions. Heavy-ion collision experiments at AGS (BNL), SPS (CERN) and RHIC 
(BNL) as well as future ones planned at FAIR (GSI), NICA (JINR), LHC (CERN) 
make it possible to explore the phase diagram of the nuclear matter in a broad pa- 
rameter range of temperature and baryon density. One of the most important issue 
addressed in nucleus-nucleus collisions is the possibility for nuclear matter to un- 
dergo phase transitions into a new state of matter. At least two of these transitions, 
chiral and deconfinement, are commonly expected. Though at zero baryon chemical 
potential both transitions seem to have the same critical temperature, it is not clear 
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at the moment whether this situation will remain at the large chemical potential 
or the the chiral phase tran sition anticipates the deconfinement one, allowing for 
the "qu arky onic" 

; tates [H2] with the massless confined quarks . Despite the r ; cent 

progress the lattice QCD calculations at finite baryon chemical potential do not 
provide fully reliable information on the subject. 

One of the most exciting signals of the deconfinement and the chiral phase tran- 
sitions in heavy-ion collisions, the chiral magnetic effect, suggested in^, predicts the 
preferential emission of charged particles along the direction of angular momentum 
in the case of the noncentral heav y- ion collisions due to the presence of nonzero 
chirality. As it was stressed in both the deconfinement and chiral phase tran- 
sitions are the essential requirements for the chiral magnetic effect to take place. 
The first one is needed, since only in this case the soft quarks can be separate d by 
a distance larger than the radius of a nucleon. The second one is also required 
since nonzero values of the chiral condensate drive asymmetry between the number 
of right- and left-handed quarks to zero. 

The effects caused by a strong magnetic field are not limited by the chiral mag- 
netic effect They include also the induced chiral symmetry breaking mod- 
ification of the nature of the chiral phase transition (e.g. turning the crossover 
pha se transition to the first-order one through influence on the c hiral cond en- 
sate ^^^^^^ J), influence on the possible color-conducting phas es j ^j^ ^ and 
the pion condensate spontaneous creation of the axial currents l -^U^ I, and for- 
mation of the 7To-domain walls'^ 1 . Recently, the effect of a large magnetic field on 
the sound velocity of a propagating plane wave was studied in Ref. 

The key quantity of these effects is a magnitude of the background magnetic 
field strength created in heavy-ion collisions. The early estimate of the magnetic 
field for the RHIC energy was made in Ref. It was shown that the field may 
reach very high values eB ~ 3 • ~ 3 x 10 18 Gausfl The aim of this article is to 
improve a qualitative estimate, make the quantitative calculation of the magnetic 
field for heavy-ion collisions at different impact parameters and different energies, 
study its characteristics. The calculations will be carried out within the microscopic 
transport model, namely the Ultrarelativistic Quantum Molecular Dynamics model 
(UrQMD), and will be complemented by analytical considerations. 



2. Magnetic field in heavy-ion collisions 

The Ultrarelativistic Quantum Molecular Dynamics model is a microscopic model 
used to simulate (ultra)relativistic heavy-ion collisons in the energy range from 
Bevalac and SIS up t o AGS , SPS and RHIC. The detailed description of the model 
can be found in Refs. E2E22 

In this paper we focus our investigation on the magnetic field evaluated in the 
center of the created fireball and on the energy density of nuclear matter in the 

a Probably the first rough estimate of the magnetic field in heavy-ion collision was made in Ref. 
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Fig. 1. The transverse plane of a noncentral heavy- ion collision. The impact parameter of the 
collision is b. The magnetic field is to be calculated at the origin O and along the y-axis. 

region surrounding this point. 

The magnetic field strength at a position x and time t is defined by the Lienard- 
Wiechert potentials 



where the fine-structure constant qem ~ 1/137, Z n is an electric charge of the nth 
particle (in units of the electron charge) and R n — x — x n is a radius vector of 
particle, v n is a particle velocity. The quantities v n and x n are taken at the time 
moment t' retarded with respect to the observation time t to be defined implicitly 
by the following equation 



where we use the natural units and set c = 1. Summation is to be carried out over all 
charged particles. However, to avoid uncertainties coming from participant contri- 
bution, the summation was performed only over the spectators in our calculations. 

In principle, the magnetic field can be calculated for any space points. We limit 
ourselves only to evaluation of the field at the central point of the reaction volume 
O (see Fig. [IJ and the dependence on y-coordinate at z = x = 0. 

On the basis of the expression (fTJ) let us draw several conclusions for the main 
properties of the magnetic field at the origin O. First of all, it is evident from 
symmetry reasons that the magnetic field will be negligible for a collision with a 
small impact parameter. From the same symmetry considerations one obtains that 
the field will have only nonzero i^-component. Second, the field will be negligible for 
low bombarding energies because the field strength is proportional to the particle 
velocity. On the other hand, for very high ultrarclativistic energies of a collision 
the contribution to the magnetic field is feasible only for particles close to the 
transverse plane (^R n — R n Vnj ~ 0. The contribution from particles away from the 
transverse plane is suppressed by the factor (1 — v 2 ). From the expression fl} it 
also follows that the characteristic magnetic field dependence on nuclei charge is 
given by eB ~ Z/R 2 , where R is the characteristic length scale proportional to the 
nuclei radius. For stable nuclei we have R ~ A 1 / 3 ~ Z 1 ! 3 and thus eB ~ Z 1 / 3 , 




(1) 



x n (t')\ +t' = t, 



(2) 
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Fig. 2. The time evolution of the magnetic field strength eB y at the central point O (see Fig[TJ in 
Au-Au collisions with impact parameter b = 4 fm in the UrQMD model, in one event ("1 ev.") and 
averaged over 100 events ("100 ev."). The symbols are plotted every At = 0.2 fm/c for E; a (, = 60A 
GeV and At = 0.01 fm/c for ^^NN = 200 GeV. 



demonstrating a weak field dependence on the nuclei charge (see similar estimate 
in Ref.QH). 

In Figs. and [31 the time evolution of the magnetic field strength for SPS 
and RHIC energies is shown. The magnetic field is created in the noncentral Au- 
Au collision with the impact parameter 6 = 4 fm. The resulting field strength is 
averaged over 100 events to reduce statistical fluctuations. It is clear, however, that 
the magnetic field in one individual event can be significant for an observable effect. 
The nonzero components of the magnetic field strength eB y in a single event are 
plotted on the same figures to demonstrate its deviations from the average value, 
see Fig. [2] In Fig. [3l the magnetic field strength eB y was estimated also assuming 
a collision of two uniformly charged Lorenz contracted noninteracting spheres with 
the radius R = 7 fm. The spheres move with the velocity defined by the collision 
energy. The results obtained in this semianalytical model and the UrQMD one agree 
to each other. The magnitude of the magnetic field estimated in Ref. ^for an earlier 
stage of Au-Au collision at the RHIC energy ^/snn = 200 GeV and the impact 
parameter b = 4 fm is about eB « 1.3 - m^, which is close to our calculations^. Note 
that this magnetic field strength is higher by about 4 orders of magnitude than that 
in the surface of magnetar It is impossible to make steady fields stronger than 
4.5 • 10 5 Gauss in the lab because the magnetic stresses of such fields exceed the 
tensile strength of terrestrial materials. 

In Fig. |4] we show the dependence of the magnetic field strength on the coordinate 
y for a fixed time corresponding to its maximal value at y — 0. As it is seen, the 
field stays approximately constant up to y ~ 5 fm, demonstrating the high degree 
of homogeneity in the central region. 



b We measure the magnetic field strength in the units of the pion mass squared, using for definitc- 
ness m-x = 140 MeV. However, the magnetic field strength can be translated into the CGS system 
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Fig. 3. The time evolution of the magnetic field strength eB y at the central point O (see Fig[TJ in 
Au-Au collisions with impact parameter, 6 = 4 fm, in the UrQMD model, for different bombarding 
energies. The symbols are plotted every At = 0.2 fm/c for E; a (, = 60A GeV and At = 0.01 fm/c 
for ^/s jviv = 200 GeV. The magnetic field obtained by modelling the gold ions as two Lorenz 
contracted non-interacting uniformly charged spheres with radius Ft = 7 fm are shown by dashed 
lines. 
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Fig. 4. The dependence of the magnetic field on the coordinate y. The calculation was carried 
out for Au-Au collisions with impact parameter 6 = 4 fm in the UrQMD model. The magnetic 
field is taken at fixed time corresponding to its maximum value at y = 0. 



The large value of the magnetic field strength itself does not guarantee possible 
observable effects, but additional requirements are needed. For the chiral magnetic 
effect the system should be in the deconfmement and the chiral restored phase. To 
demonstrate that matter in the central region is presumably in the QGP phase, we 
calculate the energy density. We define the center region of the created fireball as 
a Lorentz contracted box with the transverse and longitudinal size l x = l y = 4 fm, 
l z = (4/7 cm ) fm, where 7 cm is a Lorentz factor calculated in the center of mass 
frame. The energy density is defined as e = E/Vb ox , where Vb ox — §^/lcm fm 3 and 



by the following identity ml = 140 2 X 0.512 ■ 10 14 Gauss w 10 ls Gauss. 
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E is the energy of particles in the box 

E= £ \Jpf + m l ( 3 ) 

(:r,y,z)£box 

The calculated time evolution of the energy density in the central region is shown in 
Fig. [5] Even for the bombarding energy 2?j {, = 10 A GeV and the impact parameter 
b = 4 feu, the energy density in the central region is sufficient to reach the QGP 
phase. However, the question whether the system is in local equilibrium and can be 
described by the temperature and chemical potential(s) is still open. In the frame- 
work of the UrQMD model for central collisions the equilibration was considered in 
Refs.™. 

The magnetic field and energy density dependence on the impact parameter are 
depicted in Fig. [5] for Au-Au collisions at bombarding energies Ei a b = 60A GeV 
and y/spjN = 200 GeV . It is shown that the smaller the impact parameter b the 
larger the energy density e and the smaller the magnetic field eB y is, as it could be 
expected from simple symmetry considerations. 



3. Conclusions 

The magnitudes of the magnetic field strength play an important role in estimations 
of possible observable effects of the deconfinement and chiral phase transitions in 
heavy-ion collisions. The magnetic field strength in heavy-ion collisions at RHIC 
energies was estimated by D. Kharzeev, L. McLerran and H. Warringa following 
simple considerations, see ^ for details. For energy ^snn = 200 GeV, the authors 
obtained the magnetic field eB y m 3 • whic h is strong enough to provide nec- 
essary conditions for the effect to be measured ^Zl. In this paper we employed the 
UrQMD model for calculation of the magnetic field strength at various bombarding 
energies E lab = 10, 60, 160A GeV and ^/IWn = 130, 200 GeV. We showed the time 
evolution of the magnetic field, its dependence on the transverse coordinate y and 
the impact parameter b. The obtained results demonstrate that the fields emerging 
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Fig. 6. The magnetic field (left panel) and the energy density (right panel) dependence on the 
impact parameter b for E; a (, = 60A GeV and ^/sjvjv = 200 GeV. 

in heavy-ion collision has a high degree of homogeneity in the central region. The 
analysis of the energy density in the central region allows us to draw the conclusion 
that the considered system may reach the deconfinement and chiral restored phase. 
To all appearance, the requirements to observe the chiral magnetic effect are fully 
satisfied. We compared the UrQMD calculations with a simple semianalytical model 
considering heavy ions as two Lorenz-contracted uniformly-charged noninteracting 
spheres. The semianalytical results are in a satisfactory agreement with the dynam- 
ical ones. The former underestimate the UrQMD calculations by only few percent. 
Using this model we can estimate the lowest bound of the maximal magnetic field 
strength at the LHC energy ^/snn = 4.5 TeV to be about 15 • m% in collisions of 
Pb-Pb ions with the impact parameter 6 = 4 fm, which is twice higher than that 
estimated in Ref. 032. 
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